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ABSTRACT 

i 

The  Method  of  Auxiliary  Sources  (MAS)  is  applied  to  3D  scattering  problems  involving 
spherical  impedance  scatterers.  The  MAS  results  are  compared  with  the  reference 
spherical  wave  expansion  (SWE)  solution.  It  is  demonstrated  that  good  agreement  is 
achieved  between  the  MAS  and  SWE  results. 

INTRODUCTION 

The  Method  of  Auxiliary  Sources  (MAS)  is  a numerical  technique  applicable  to 
electromagnetic  scattering  problems.  In  the  general  case,  a set  of  spatially  impulsive 
electric  and/or  magnetic  sources  is  introduced  to  radiate  an  approximation  to  the 
unknown  scattered  field.  These  so-called  auxiliary  sources  are  located  on  an  auxiliary 
surface,  typically  conformal  to,  and  enclosed  within,  the  physical  surface  of  the 
scatterer.  Point  matching  of  the  boundary  condition  on  the  physical  surface  is  enforced 
to  determine  the  complex  amplitudes  of  the  auxiliary  sources.  The  MAS  originates  from 
an  application  of  a special  case  of  the  Method  of  Moments  (MoM),  utilising  spatially 
impulsive  expansion  and  testing  functions,  to  a generalised  surface  integral  equation 
formulation  [1],  An  overview  of  MAS  is  given  in  [2],  Utilisation  of  MAS  for  numerical 
solution  of  various  2D  scattering  problems  and  3D  PEC  and  dielectric  scattering 
problems  has  been  reported  earlier  [3],  [4],  [5],  The  purpose  of  this  work  is  to 
investigate  the  performance  of  MAS  when  the  method  is  applied  to  scattering  problems 
involving  spherical  impedance  scatterers.  The  spherical  wave  expansion  (SWE)  solution 
is  developed  and  used  as  a reference.  The  scattering  problem  under  investigation  is 

illustrated  in  Figure  1,  together  with  the 
introduced  Cartesian  and  spherical  co-ordinate 
systems.  The  spherical  scatterer  of  radius  ro  is 
illuminated  by  an  x-polarised  uniform  plane  wave 
(E',H')  of  wavelength  X and  propagating  in  the  z- 
direction.  The  Standard  Impedance  Boundary 
Condition  (SIBC)  holds  on  the  physical  surface  of 
the  scatterer,  rx(Exf)  = £rxH.  The  chosen 
auxiliary  sources  are  pairs  of  crossed  0 - and  <p  - 
directed  electric  Hertzian  dipoles  located  on  a 
sphere  of  radius  r0-S.  The  auxiliary  source  pairs, 
as  well  as  the  matching  points,  are  placed 
equidistantly  in  the  angular  co-ordinates  ( 6,<p ). 
The  total  number  of  source  pairs  is  denoted  by  N. 


Figure  1 : The  scattering  problem  geometry 
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The  SWE  solution  is  given  by  (see  [6,  Chapter  2]  for  definitions  of  symbols) 
E = kj \jrj  X Qs  ,m,ii^s,m,n  »H  =-ikJrj  X Qs,m,n^^s  mn’ Qs'Jiim  = & a/b. 


s,m,n 


a = in+] -Jn (In  + 1)(-1)4  5 m R(^}s  nR(^n{-8 s]{8 nuX  +8m-\)  + 8s,2^m-\  ~^m,\))~ 

in+2^K(2n  + l)(-8s^yVn  (SmA  + 8mt_0  + - SmA )), 

b=-k/jj  (-1)3-'"--9  + oyn-if1  (R™sy , 

*g  BhUHkro),  RVy/n(kr0),  R™mh2hkn>). 


NUMERICAL  RESULTS 

We  define  the  boundary  condition  error  (BCE)  by 

I M,N 

BCE  ^ x I *m,n  x (E,„,„  x r ) - £ r x H„,  I (1) 

MN  IE'  I m,n 

The  summation  is  performed  over  the  scatterer  physical  surface.  Figure  2 shows  the 
BCE  as  a function  of  N for  four  different  scatterer  surface  impedances  £ and  three 
values  of  £(♦:  8/X=0.2, 1 : 8/X=0.5,  s : S/X=0.8).  In  all  cases,  r0l~k-\  is  chosen.  It  is 
observed  that  the  BCE  attained  for  8/1=0. 5 is  the  lowest,  and  that  it  starts  with  a rapid 
decrease,  whereafter  it  attains  a constant  level. 


Figure  2 
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In  Figure  3,  the  normalised  bistatic  radar  cross  section  (BRCS)  results  obtained  using 
MAS  for  four  different  scatterers  and  for  two  different  planes  of  observation  (<p=0  and 
qp=7i/ 2)  are  compared  with  the  corresponding  SWE  results.  In  all  cases,  ^=100+jl00L>  is 
chosen.  There  is  a very  good  correspondence  between  the  MAS  and  SWE  results.  The 
small  discrepancies  can  be  diminished  further  if  N is  increased. 


Figure  3 


CONCLUSION 

It  is  seen  that  the  initial  decrease  in  the  boundary  condition  error  of  the  MAS  numerical 
solution  is  relatively  fast.  Also,  the  limiting  error  level  achieved  is  small  enough  for  the 
MAS  to  be  able  to  produce  good  approximations  to  exact  results  at  a relatively  small 
computational  cost.  However,  the  existence  of  a finite  limit  level  for  the  MAS  boundary 
condition  error  indicates  that  the  present  numerical  implementation  of  the  method  is  not 
convergent.  Bistatic  radar  cross  section  results  obtained  by  MAS  are  found  to  be  in 
good  agreement  with  the  reference  spherical  wave  expansion  (SWE)  solution. 
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